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We describe a random matrix model suitable for the sim- 
ulation of the eigenvalues of the Dirac operator on the lattice 
for Wilson fermions. We compare the obtained global eigen- 
value spectrum for various values of the hopping parameter k 
with lattice results of Kalkreuter. The agreement is surpris- 
ingly good. 
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Recently, it became clear that the microscopic spec- 
tral properties of the lattice QCD Dirac operator are 
universal and can be reproduced by simple models that 
only share basic symmetry properties with real QCD 
[|l|,||. Such models are provided by random matrix the- 
ory (RMT) This universality has recently been 
demonstrated for the staggered lattice Dirac operator 
in quenched ^ and unquenched [|| SU(2). Since the 
Banks-Casher formula S = 7rp(0)/V ||] links the spec- 
tral density at zero virtuality to the chiral condensate E 
the distribution of the small eigenvalues is of great im- 
portance for, e.g., the understanding of the chiral phase 
transition. RMT has also solved a long standing problem 
of lattice calculations at finite chemical potential |^ . Re- 
cently, also the predictions for the energy scale at which 
RMT and lattice QCD start to deviate was confirmed, 
which among others tested the validity of the Gell-Mann- 
Oakes-Renner relation on the lattice ||]. While RMT 
makes reliable predictions only for microscopic spectral 
fluctuations, these successes encourage us to see if the 
global spectral properties can be reproduced. This is the 
aim of our contribution for the case of Wilson fermions. 
Since the universality argument only applies to micro- 
scopic properties one cannot hope to fit global spectra 
except for very special cases. We analyze such a case, 
namely, gauge theories with infinitely strong coupling, 
i.e., P ^ 0. Such systems can be hoped to be sufficiently 
chaotic as to show random matrix characteristics even 
on global scales. Furthermore our model studies suggest 
that to describe Wilson fermions a 4 x 4 block structure is 
needed in random matrix theory. We expect this to hold 
true also for a description of microscopic properties at 
non-zero /3. Random matrix models for Wilson fermions 
at non-zero P do not yet exist, but are of great practical 
importance. 

We start our analysis with the Euclidean action of lat- 
tice SU(2) theory with Wilson fermions in the fundamen- 
tal representation which can be written as 
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with the Wilson parameter r that we set to 1 in the fol- 
lowing [^. The gauge fields are contained in the link 
variables U. With this action a gauge invariant partition 
function can be constructed from which one can obtain 
vacuum expectation values of operator products in the 
usual way. In these partition functions we average over 
all gauge field configurations. 

In random matrix theory we substitute the Dirac op- 
erator which includes the gauge fields by random ma- 
trices of a particular ensemble to model the very strong 
fluctuations of the Dirac operator when calculated with 
lattice gauge theory. The integration is then performed 
over the independent entries of the matrices [|l|. With 
this approximation the gluons decouple from the quarks 
and can be integrated out, i.e., they can be neglected in 
random matrix models if one is only interested in quark 
observables. The symmetry properties of the random 
matrix depend on the underlying gauge group and the 
fermion representation. Usually in QCD we are dealing 
with fermions in the fundamental representation of the 
SU(3) gauge group in which case one has to use the Gaus- 
sian Unitary Ensemble (CUE). In order to compare our 
results with Kalkreuter who investigated the opera- 
tor 75 + m) for massive Wilson fermions in an SU(2) 
gauge field background we have to use matrices of the 
Gaussian Orthogonal Ensemble (GOE) 0. The Eu- 
clidean partition function with two random matrices can 
be written according to the above arguments as 



V[A,B]e- 
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with the parameters T^a and that scale the distribu- 
tion variance of the Gaussian Ensembles. We will now 
specify the operator ^ + m in Eq. (||) and explain why 
we are using two different random matrices. 

In the following we separate the Dirac spinors into 
left- and right-handed fields V^l,r = ^(1 T Ts)'^ with 
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the Euchdean 75 = —7^. Furthermore we distin- 
guish between even and odd lattice sites. These 
two sub-lattices couple rather independently (the two 
groups are often called 'red' and 'black' lattice sites 
in analogy to the colors of a checkerboard). The 
spinor field is then written as (V^r, V^l? V^r? V^l)' where 
'^R°L ~ (V^R^ (-'-)' • • • ' V^r^It)) vectors with respect 
to the lattice sites . . . , x^{N). 

In the above basis the mass term of (|l|) is simply a 4 x 4 
block diagonal matrix and the interaction term becomes 



-^(V'£V£+V'aVji + ^£V£ + ^JiVa) (3a) 

(3b) 
(3c) 

- ^5«(-^lV^r - ^rVC - ^lMVr - ^^^MVD- 

(3d) 

The notation /i indicates that /i is not a free Lorentz 
index but is contracted with a corresponding one that is 
hidden in the ?/^'s. We have expanded the link variables 
to first order in A^, i.e., [/^ = 1 + i^aA^. The term 
then leads to the constant block matrix 
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In ( p^ ) we replace with i/ ga times the random matrix 
A such that we get the block matrix 
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(4b) 



In (^) we do the same with ^ but calling the random 
matrix B which gives us 
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In the remainins; term (||) which has the same structure 
as ( |3d| ) we are faced with a current-like form, bilinear 
in the fields. As in (|^) different 7 matrices are used 
depending on the lattice sites which makes this term suf- 
ficiently random to be absorbed in the matrix B. To ob- 
tain the above matrices we have rescaled the Dirac fields 
by 



The expansion ofUfi we used to argue for our assumed 
block structure is justified in the perturbative domain 
of small g. Here we are rather interested in the ergodic 
regime of fully developed chaoticity, which corresponds 
to rather large g. However, we are only interested in the 
general symmetry properties, which are the same in both 
regimes. 

Adding up the above contributions (||) we arrive at the 
following expression for the complete Dirac operator: 



( 1/2K I- A B \ 
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(5) 



Neglecting the distinction between even and odd fields 
and looking at the chiral structure only. If) has the struc- 
ture of a 2 X 2 block matrix, namely, ( 2 ) since it is 
anti-commuting with 75 . This ensures that the eigenval- 
ues of Jf) are distributed around zero symmetrically. Let 
us stress that we have in no way derived Eq. (P). We 
basically only give some hand-waving arguments for the 
form in Eq. (^) which is then tested against the numeri- 
cal data. Definitely in the end more realistic models will 
differ from (||) but we believe that they will share the 
4x4 structure. 

Kalkreuter investigated in his analyses the operator 
75 (^ + m) so we have to multiply (|^) by 75. Here we 
must carefully keep in mind that the B matrix resulted 
from the ^ term which implicitly contains the 7 matri- 
ces. When we also define our spinor basis according to 
('^R' ~'^L5 '^R' '^l) finally arrive at 



75(^ + m) = 
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At first we tried to determine the spectral density of 
the above operator analytically. For this purpose we 
chose (|^) as the matrix in the determinant of (||) and 
integrated out the 2 • ^y(y + 1) real matrix variables 
by means of the Hubbard-Stratonovich transformation 
which lead to a partition function with integration over 
only 8 complex variables. It turned out that the resulting 
saddle point equations are too complicated to be solved 
analytically so that we decided to diagonalize the opera- 
tor matrix numerically. 

We calculated the eigenvalues of this matrix for a large 
number of random variables matching the selected en- 
semble (GOE). We display the results in histograms for 
several values of the parameter n = (2m + 8)~^ (see 
Fig. [^). Our results are on the left side of this figure. 
They are compared with the corresponding /3 = data 
of Kalkreuter taken from [|l2| which are shown on the 
right side. 
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FIG. 1. spectral densities of the massive Dirac operator for various values of the hopping parameter k. On the right hand 
side are shown Kalkreuter's results for /3 = 0.0 extracted from the first picture of Fig. ^ but with eigenvalue normalization as 
on the left hand side. The chosen parameters on the left are = 2/25, — 8/25. The hopping parameter is k — 1/2, 1/4, 
1/5, and 1/8, respectively. 
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One can clearly see the splitting of the spectral density 
p{X) into two distinct symmetric parts at a critical value 
of hi between 0.2 and 0.25 as well as a slight dent on 
top of the half-circle like densities. The structures come 
about as follows. The separation of the two half-circular 
structures is due to the diagonal term 1/2K,. The splitting 
of the two peaks for each half-circle like structure results 
from the constant '1' in (||). These dependences are in 
good agreement with the data. The matrices A and B 
fluctuate according to the random matrix constants T^a 
and (see (|^)). If T^a and T>b are large the structure 
gets washed out. If and are different the half- 
circle like structures become asymmetric. 

Since the gauge fields in the action are completely re- 
placed by random variables the best agreement has to be 
expected for /3 = 0. The comparison shows that in this 
limit of extremely strong coupling RMT is actually able 
to fit the global eigenvalue spectrum surprisingly well. A 
detail which is missed is the slight difference in height of 
the two maxima for, e.g., positive A. Also the falloff for 
large eigenvalues is steeper for the lattice results. 

With decreasing coupling strength the agreement for 
the global spectrum becomes worse, as is seen by compar- 
ing the left hand side of Fig. |l| to the lower two pictures 
of Fig. ^, and one reaches the usual situation, in which 
RMT can only describe the microscopic fluctuations. 

Let us conclude: We defined for the first time a ran- 
dom matrix model which is suitable for the description 
of the eigenvalue spectrum of the Dirac operator (mul- 
tiplied by 75) for an SU(2) gauge theory with Wilson 
fermions. A similar model should allow to study, e.g., 
the distribution properties of the lowest eigenvalues for 
SU(3)-Wilson fermions and for the operator (p-{-m) 
instead of 75 (^ + m), which is a problem of great prac- 
tical importance. 
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FIG. 2. Kalkreuter's results for /3 = 0.0, 1.8, and 2.8, re- 
spectively. The eigenvalues are normalized to lie between —1 
and 1. The symbols +, x, O, □, and A correspond to 
quenched data with k values of 0.25, 0.20, 1/6, 0.15, 0.125, 
and 0.10, respectively. The solid line shows unquenched data 
at At = 0.15 and (3 2.12. 
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